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*Pht«  report  describes  an  analysis  of  «  model  of  the  active 
sonar  detection  process.  The  analysis  was  done  In  order  to 
estimate  the  effects  of  signal  excess  fluctuations  on  the 
probability  of  detecting  a  sonar  echo.  The  analysis  la  an 
extension  of  an  analysis  that  Is  described  in  Reference  l. 

The  model  has  the  following  characteristics:  Each  retxirn 
corresponding  to  a  resolution  cell  has  an  associated  signal 
excess.  A  return  la  condltlcmally  classified  as  an  echo  If  and 
only  If  the  signal  excess  essoclated  with  the  rettim  Is  greater 
than  or  equal  to  sero  (a  success).  A  return  la  unconditionally 
classified  as  an  echo  If  and  only  If  It  la  the  success  In  a 

sequence  of  no  greater  than  n  returns.  This  detection 
criterion  has  been  referred  to  as  a  k^*out“Of-*n  criterion. 

For  a  resolution  cell  that  does  not  contain  a  target,  the 
signal  excess  Is  n^ative  infinity.  F.*>r  a  resolution  cell  that 
does  contain  a  target,  the  signal  excess  is  a  rand(^  rjuantlty 
that  Is  determined  by  a  stochastic  process.  Three  stochastic 
processes  were  considered  in  the  analysis.  The  first  is  the 
gauss'-markov  process^  the  second  is  the  laabda^-slgma  jump  process 
and  the  third  is  a  acdiflcatlon  of  this  process.  In  the  analysis 
that  is  described  in  Reference  1,  only  the  lai^a-sigi&a  jimp 
process  ajid  its  modificatlou  were  considered.  A  discussion  of 
the  gauss -Barkov  process  which  supports  its  inclusion  in  the 
analysis  is  contained  in  Reference  2  and  an  example  of  its  use  in 
an  active  sonar  detection  model  is  contained  in  Reference  3. 
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TbM  active  sonar  detection  models  that  are  described  in 
Beference  3,  E^ference  4  and  Eeference  5  are  all  based  on  a 
h*^iit-*of-n  criterion.  In  Eeference  S,  it  is  suggested  that  k 
«3ual  to  three  and  n  equal  to  five  give  an  adequate  description 
of  the  active  sonar  detection  process  in  vhich  an  operator 
determines  the  classification  of  the  sonar  returns,  in  Eeference 
6,  arguBsntit  based  on  ^lerational  data  are  given  both  for 
determining  sonar  return  classification  in  terms  of  signal  ercess 
and  for  a  k-out'>of~n  criterion. 

Since  each  return  correepondii^  to  a  resolution  cell  has  an 
associated  signal  excess,  the  model  can  accommodate  false 
contacts  cased  on  signal  inputs  but  not  falsa  contacts  based  on  . 
noise  alone.  However,  if  in  the  model  the  probability  of 
conditionally  classifying  noise  as  an  echo  were  .001,  then  with  a 
3“Out-of-5  criterion  the  chance  of  t  false  contact  based  on 
noise  alone  would  be  less  than  one  in  ten  million. 


II.  ‘SSbM  kstivm  Sonar  Dataction  Model 

In  the  aodei  of  the  active  sonar  detection  process,  the 
claseification  of  an  echo  is  determined  by  its  signal  esccess.  A 
sonar  return  is  conditionally  classified  as  im  echo  if  emd  only 
if  the  signal  excess  is  greater  than  or  equal  to  zero.  The 
signal  excess  is  the  difference  between  the  echo  signal ~to~noise 
ratio  in  decibels  and  the  sonar  receiver  recognition 
dif f erentizU. ,  both  of  which  are  independent  random  variedsles  in 
the  model.  With  Xse(^)  signal  excess  for  an  echo  generated 

at  a  time  indicated  by  the  index  t,  the  signal  excess  is: 

(1)  XsE(t)  -  10  log(S/N)  -  XRD(t) 

where  10  log(S/H}  is  the  signal -to-noiee  ratio  in  decibels  for 
an  echo  generated  at  that  time  and  X|^(t)  is  the  recognition 
differential  for  the  echo.  The  random  signal-to-noise  ratio  in 
decibels  is  determined  by  the  following  sonar  equation: 

(2)  10  log(S/N)  -  Xsi,(t)  -  2  X^LCt)  -f  XTS(t)  -  tX^LCt)  -  Xoi(t)] 
where,  Xg^Ct)  is  the  sonar  source  level,  X«j»L(t)  is  the  one¬ 
way  transmission  loss  between  the  target  and  the  sonar,  X^s(t) 
is  the  target  strength,  XjjeC^)  is  the  noise  level  and  XQj(t) 
is  the  sonar  directivity  index.  The  random  variable  Xjji,(t) 
accounts  for  ambient  noise,  self  noise  and  reverberation  and 
represents  a  power  sum  of  these  quantities. 

From  Equations  1  and  2,  the  mean  value  of  the  signal  excess 
is  given  by: 

(3)  SE(t)  -  SI,(t)  -  2  TL(t)  TS(t)  -  tNL{t)  -  DI(t)]  -  RD(t) 
Where  each  term  represents  the  mean  value  of  its  corresponding 


rtsMoa  variable.  In  the  model,  the  mean  values  are  the  values 
that  vould  be  used  in  the  sonar  equation  that  Elation  3  usually 
represents,  from  Equation  l  and  Equation  3  the  signal  excess  can 
be  e^^essed  as  follows: 

(4)  Xse(^)  “  SE(t)  +  X(t) 

where  X(t)  is  a  random  variable  whose  mean  is  zero  and  whose 
varijmce  is  equal  to  the  variance  of  XsE(t).  Equation  4  implies 
that  for  each  index  t  that  is  involved  in  describing  an 
encounter  there  is  a  raxuiom  variable  X(t)  that  determines  the 
fluctuations  in  the  signal  excess  at  the  time  corresponding  to 
the  index  t. 

The  three  stochastic  processes  that  were  used  in  the 
analysis  to  determine  the  random  variables  represented  by  X(t) 
are  described  in  Section  III.  Unfortunately,  it  appears  that  the 
following  quotation  from  Reference  7  has  seme  relevance  for  each 
of  the  processes:  "As  discussed  in  Volume  I,  it  is  common 
Practice  Ut  performance  analyses  to  model  signal  or  noise  or 
signal  excess  fluctuations  as  stochastic  processes.  The  choice 
of  a  specific  process  may  be  based  on  experience  but  usually 
tends  toward  mathematical  convenience." 
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in.  Vvo  Siaalation  Prograas 

Tvo  simulation  programs  for  an  IBH-PC  were  used  in  the 
analysis  of  the  active  sonar  detection  model.  The  two  programs 
differ  in  the  stochastic  process  that  determines  the  reuidom 
components  of  the  signal  excess.  In  the  first  program,  the 
stochastic  p^’ocess  is  either  a  leuabda-sigma  jump  process  or  a 
modification  of  a  lambda^-sigma  jump  process.  In  the  second 
program,  the  stochastic  process  is  a  gauss-markov  process.  The 
progreuos  which  are  written  in  BASICA  are  described  in  Appendix  1 
euid  are  listed  In  Appendix  2. 

The  programs  simulate  encounters  between  a  submarine  and  an 
active  sonar.  In  an  encounter,  the  sonar  source  level,  the 
tremsmlsslon  loss,  the  submarine  target  strength,  the  noise  level 
and  the  «onar  directivity  index  are  independent  of  azimuth.  In 
addition,  the  course,  speed  and  depth  of  the  submarine  and  the 
depth  of  the  transmitting  and  receiving  hydrophones  of  the  sonar 
are  constant  throughout  an  encounter.  This  implies  that  the 
track  of  the  submarine  relative  to  the  sonar  is  a  straight  line 
and  that  the  encounter  can  be  determined  by  specifying  following 
quantities:  the  depth  of  the  submarine,  the  depth  of  the 
transmitting  hydrophone,  the  depth  of  the  receiving  hydrophone, 
the  relative  speed  of  the  submarine,  the  horizontal  range  of  the 
submarine  at  its  closest  point  of  approach  (CPA) ,  the  time  at  the 
beginning  of  the  encounter  and  the  time  at  the  end  of  the 
encounter.  An  encounter  of  this  kind  that  begins  and  ends  at 
ranges  beyond  which  the  probability  of  detecting  the  submarine  is 
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essentially  zero  is  called  a  complete  straight  line  encounter. 

?or  a  ooB^lete  straight  line  encounter,  only  the  horizontal  range 
at  CPA  is  required  to  determine  the  encounter  geometry,  ^^lis 
range  is  called  the  lateral  range.  The  probability  that 
detection  occurs  during  a  complete  straight  line  encounter  can  be 
emresaed  as  a  function  of  the  lateral  range.  The  function  is  ' 
called  the  lateral  range  function  or  lateral  reuige  curve. 

The  lateral  range  functions  plots  that  are  in  Section  IV 
were  generated  with  the  two  simulation  programs.  The  mean 
transmission  loss  values  that  were  used  to  this  are  plotted  in 
Figure  1  of  Section  IV.  The  values  correspond  to  horizontal 
remges  from  2000  meters  to  100,000  meters  in  2000  meter 
increments.  For  a  target  range  that  is  not  an  integer  multiple 
of  2000  meters,  the  me2ui  transmission  loss  value  that  is  used  in 
the  programs  is  the  value  corresponding  to  the  greatest  integer 
multiple  that  is  lees  than  the  target  range.  Relative  to  the 
programs,  the  transmission  loss  is  determined  by  a  step  function 
with  a  step  every  2000  maters. 

Because  of  the  step  function  nature  of  the  transmission 
loss,  the  relative  track  of  the  target  in  an  encounter  is  divided 
into  sectors  such  that  when  the  target  ie  within  a  sector  the 
transmission  loss  is  constant.  The  simulation  programs  generate 
estimates  of  both  the  conditional  and  the  unconditional 
probability  that  detection  occurs  on  a  sector. 

Clearly,  the  encounters  represented  by  the  simulation 
programs  ar©  abstractions.  In  particular,  except  for  the  mean 
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t3afcni»i«iicn  loaNS  TL(t) ,  all  tenas  on  the  right  side  of 
Hguaticm  2  asm  constant  during  an  encounter.  Changes  due  to 
orisntation,  location  and  relative  motion  are  ignored.  For 
example,  in  addition  to  ignoring  the  change  in  the  noise  field 
with  azimuth,  the  change  in  the  noise  field  due  to  the  doppler 
gain  that  resi’  ts  from  the  radial  notion  of  the  sonar  and  the 
eubiBarine  is  also  ignored.  (More  detailed  models  are  described 
in  Reference  3,  Reference  4  and  Reference  5.)  However,  the 
simulation  prograiss  were  assumed  to  be  adequate  to  satisfy  the 
purpcee  of  the  analysis. 

Because  of  the  conditions  described  above.  Equation  3  can  be 
written  as: 

(5)  SE(t)  -  SL  -  2  Th(r(t)]  +  TS  -  (IJL  -  DI)  •  RD 
where  the  transmission  loss  at  a  time  corresponding  to  the  index 
t  Is  determined  by  the  sonar  transmitting  hydrophone  depth,  tdse 
sonar  receiving  hydrophone  depth,  the  submarine  depth  and  r(t), 
the  submarine  horizontal  range  at  that  time. 

Both  the  lambda-sigma  jump  process  and  the  gauss-aarkov 
process  have  been  used  to  predict  operational  performance  of 
sonar  systems.  This  is  not  the  case  for  the  modified  lambda- 
sigma  jump  process.  The  modification  was  introduced  in  order  to 
dmal  with  the  observations  that  are  reported  in  Reference  8.  A 
modification  based  on  a  log-normal  distribution  rather  than  a 
shifted  rayleigh  distribution  would  have  accomplished  this  also. 
The  choice  of  a  rayleigh  distribution  was  based  on  mathematical 
convenience . 


For  a  raylalgb  randoa  variable  Y,  tba  density  function  of 
its  diatributiott  las  fy(y}  2*o*yoG3£p(‘^*y*)  wlisre  y  ^  0. 

0  *■  r(3/2),  tha  iscsn  of  tha  distribution  is:  «»  $/ak 

and  tStm  variance  la:  »  (1  -  ^*)/a. 

TbQ  isodified  laabda<>aigma  jus$>  procaaa  is  defined  as 
follova:  If  X(t)  is  a  random  variable  that  is  datensinod  by  a 
sodified  laubda-eigaui  jump  procase,  than  X(t)  «•  ~  %{t) 

vhara  H(t)  is  determined  by  a  Poisson  process  vith  a  s@an  ra-^ 
laabda  aiul  the  random  variables  Yx,  Y^,  *•*  are  independent 
rayleigh  random  veuriables  each  with  parameter  a.  I^is  implies 
that  ux  ■  0  Ox  *  *^Y«  a  »  (1  -  p®)*(l/o*),  then 

Ox  **  o.  Conaequcuitly,  if  the  randena  component  of  the  signal 
excess  is  determined  by  a  modified  lambda~6igx:a  jump  process, 
than  the  standax^d  deviation  of  the  signal  excess  can  be  made 
equal  to  sigma  by  requiring  that  a  -  <1  -  0®)»{Vo'*)« 

Iha  correlation  coefficient  between  two  random  variables  is 
a  measure  of  their  dependence.  For  the  random  variables  X{tx) 
and  K(tj}  which  represent  the  random  cca^nent  cf  the  signal 
excess  at  timaa  and  tj  and  for  tj  >  t^,  tha  correlation 

coefficiont  has  the  form  exp[-{t-^  -  tx)/T]  for  the  gauss-markov 
process,  the  lambda''sigQta  jump  process  and  the  modified  lambda- 
sigma  jta^  process.  For  each  process,  a  value  for  the  parameter 
T  and  a  value  for  the  parameter  a  ara  sufficient  to  define  the 
process.  Tho  pareineter  t  is  referred  to  ss  tlio  relaxation 
time.  It  is  the  reciprocal  of  the  parameter  lambda  in  tha 
l&Ebda-sigssia  jump  processes.  For  booh  the  laabda-sigma  jump 


jmp  processes.  For  both  it  and  the  modified  lambda-sigma  jump 
process,  t  is  the  expected  time  between  jumps.  If  t  »  0  for 
any  one  of  the  processes,  the  random  components  of  the  signal 
excess  are  independent.  Otherwise,  the  random  components  are 
dependent.  If  (1/t)  =*  0  for  any  of  the  processes,  the  random 
components  are  completely  dependent.  ‘This  implies  that  the 
random  components  of  the  signal  excess  associated  with  an 
encounter  all  have  the  same  value.  Consequently,  if  the  value  of 
one  of  the  random  components  of  the  signal  excess  is  determined 
at  some  point  in  the  encounter,  the  value  of  all  of  the  other 
random  components  are  also  determined.  Reference  9  contains  a 
discussion  of  factors  affecting  relaxation  times. 

An  acceptable  choice  for  a,  the  standard  deviation  of  the 
signal  excess,  appears  to  be  more  easily  agreed  upon  than  an 
appropriate  choice  for  t,  the  relaxation  time.  However,  the 
following  quotation  from  Reference  10  is  relevant  both  to  these 
choices  and  to  the  choice  of  an  appropriate  stochastic  process: 
"As  of  this  date  [1974]  there  is  no  unanimity  of  opinion  as  to 
which  process  [lambda-sigma  jump  or  gauss-markov]  to  use  in  a 
detection  model  or  what  values  of  the  parameters  lambda  and  sigma 
to  use.  It  is  not  even  clear  (certainly  not  to  everyone)  that 
this  kind  of  detection  model  is  uniformly  valid.  Some  half¬ 
hearted  attempts  at  ‘validating’  this  kind  of  model  have  been 
made  but  the  results  can  bo  characterized  as  at  best  fragmentary 
and  inconclusive.  In  my  opinion,  what  is  needed  is  a  full-blown 
exercise  whose  sole  purpose  is  to  investigate  the  validity  of 
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various  hypothesized  sodels  for  detection.  Past  attempts  at  this 
have  merely  been  'esc  post  facto*  eadmaation^i  of  some  sketchy 
exercise  data.** 
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2V«  liJjBiilatiosi  Seisalts 

Bnc»unter  stssilation  data  are  presented  in  this  section. 

^e  data  i^icth  are  presented  as  lateral  range  function  plots  and 
as  range  sector  probability  of  detection  tables  are  from 
simulations  that  are  based  on  the  propagation  loss  values  that 
ore  plotted  in  Fi.gure  1,  the  encounter  values  that  are  listed  in 
Table  1  and  a  3-out--of -5  criterion. 

For  the  siikUlated  encounters,  the  lateral  range  function 
plots  in  Figure  2  through  Figure  6  Indicate  the  effect  on  the 
encounter  probability  of  detection  of  the  transmission  loss 
function,  especially  the  increase  witii  range  because  of  the 
convergence  zone  component.  3ut,  in  particular,  the  figures 
indicate  the  effect  on  the  encounter  probability  of  detection  of 
a  gauss-markov  process  relative  to  that  a  lambda-sigma  jump 
process  for  various  relar:ation  times.  Figure  7  is  representative 
of  tile  affect  on  encounter  detection  probability  of  a  modified 
lambda-sigma  jump  process  relative  to  that  of  the  unmodified 
process.  For  additional  comparison  plots  for  tJie  sane  simulated 
encounters,  see  Reference  l. 

Except  for  the  CPA  range  sector,  each  range  sector  is 
associated  with  two  track  segments  that  are  symmetric  about  the 
CPA.  The  CPA  range  sector  is  associated  with  a  single  track 
segment  that  is  bisected  by  the  CPA.  Table  2  and  Table  3  lisr 
the  probability  of  detection  estimates  cn  a  track  segment  as  a 
function  of  the  track  segment's  distance  from  CPA  and  the  range 


aMoci&ted  y^ith  th«  track  segatant.  Tabla  4  and  Table  5 
dc  tide  for  l^e  conditional  probability  of  detection  estimates, 
^o  values  in  Table  2  and  Table  3  are  estimates  of  the 
probability  that  detection  will  occur  on  a  track  segment.  The 
values  in  Table  4  and  in  Table  5  are  estimates  of  the  probability 
that  detection  will  occur  on  a  track  segment  given  detection  has 
not  already  oociirred  on  an  earlier  track  segment. 

The  length  of  the  simulated  straight  line  encotmters  that 
were  used  in  the  analysis  was  restricted  by  the  requirement  that 
the  rangi  from  the  sonar  at  the  end  points  of  an  encoinnter  be 
less  than  or  equal  to  lOO  kilometers.  Because  of  this 
restriction,  the  encounter  track  length  is  120  kilometers  for  an 
encounter  lateral  range  of  30  kilometers  and  the  encc‘jat©r  track 
length  is  appromimately  200  kilometers  for  an  encounter  lateral 
range  of  10  kilometers.  Howe’»er,  for  the  simulated  acoustic 
conditions,  the  transmission  lose  for  a  range  of  100  kilometers 
is  94.2  dE.  If  transmission  Icjs  values  for  ranges  greater  than 
100  kilometers  were  ccnstant  and  aqual  to  94.2  dS  rather  than 
generally  increasing,  for  an  ©"bo  from  c  rarget  that  was  at  a 
range  greater  than  100  kilomic^rs,  the  piobability  of 
conditionally  classifying  the  echo  as  an  echo  would  be  ii.  the 
neighborhood  of  .000015.  Inus,  increasing  the  length  of  the 
simulated  encounters  should  only  sr'ik©  a  negligible  increase  in 
the  value  of  the  encounter  probability  of  detection  estimates. 

'ibe  probability  estimates  genarated  by  the  programs 
represent  samples  from  binomial  distributions  to  the  degree  of 


X2 


randcameos  of  tho  program  siaulatlona.  Howavor,  the  saaple  else 
of  th%  estlaataa  is  not  constant,  in  a  sinulatsd  oncountsr,  the 
encounter  is  terminated  %dien  a  detection  occurs.  Conse(suently, 
for  a  given  number  of  repetitions  of  an  encounter  simulation,  the 
sample  size  for  a  range  sector  depends  on  the  number  of 
detections  that  have  occurred  in  track  segments  prior  to  the 
track  segment  associated  with  the  range. sector.  Because  of  this, 
each  simulated  encounter  vas  repeated  1000  times  in  order  to 
generate  the  lateral  range  function  plot  emd  sector  probability 
of  detection  estimates.  This  number  of  samples  is  statistically 
sufficient  for  the  purpose  of  the  analysis  reported  here. 
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Figure  1.  A  plot  of  the  encotmter  transaission  loss  valvtos  In  dB 
against  their  corresponding  ranges  in  kilometers.  The  values  are 
listed  in  Appendix  4. 
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Rttccgnltion  Differential 

20 

dB 

Source  Level 

215 

dB 

Target  Strength 

5 

dB 

Hoiee  Level 

65 

dB 

Directivity  Index 

20 

dB 

Standard  Deviation:  a 

8 

dB 

Ping  Cycle  Time 

61 

seconds 

Table  1.  The  parameter  values 

for  the  simulated 

encounters 

Siliematew 


Q  4-  Hi 

Figure  2.  Lutarfil  range  function  plot*  for  a  gau»s>'Qarkov 
process  and  a  Xaiabda-sigsa  jusp  process  for  a  relaxation  time 
t  of  3  minutes. 
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Figure  3.  Lateral  rai^«  function  plots  for  a  gauss-xarkov 
procass  and  a  lafilxia~aigaa  jusp  procass  for  a  relaxation  time 
T  of  10  xinutaa. 
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Figure  5.  Lateral  range  function  plots  for  a  gauss-markov 
process  and  a  laobda-sigaa  juQp  process  for  a  relaxation  time 


T  of  zero. 


a 


Figure  6.  liateral  range  function  plots  for  a  gauss-aarkov 
process  and  a  lambda-sigma  jinap  process  for  a  relaxation  time 
T  of  infinity. 
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Figure  7.  Lateral  rang©  function  plots  for  a  leuabda-sigma  juEp 
process  and  a  modified  iaaixia-sigisa  jump  process  for  a  relaxation 
time  T  of  10  minutes. 


Table  2.  Range  actor  probability  of  detection  estimatee  for  a 
gauss-sarkov  process,  a  lateral  range  of  10,000  asters,  a 
relaxation  tine  t  of  3  ninutes  and  a  3-oot“Of-5  criterion. 
Range  sector  and  distance  froa  CP?*  values  are  in  aeters  and  zero 
probability  values  are  not  listed. 
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probability 

range  sector 

distan 

.01 

76000 

-74330 

.132 

74000 

-72312 

.089 

72000 

-70292 

.062 

70000 

-68272 

.029 

68000 

-66250 

.045 

66000 

-64226 

.007 

64000 

-62201 

.003 

62000 

-60175 

.001 

50000 

-47969 

.001 

480QG 

-45924 

.001 

40000 

-37697 

.002 

38000 

-35623 

.001 

26000 

-22913 

.007 

12000 

-4583 

.007 

10000 

0 

.001 

12000 

4583 

.001 

32000 

29343 

.001 

38000 

35623 

.001 

46000 

43875 

.002 

50000 

47969 

.001 

52000 

50010 

.001 

58000 

56116 

.002 

60000 

58146 

.007 

62000 

60175 

.040 

64000 

62201 

.030 

66000 

64226 

.046 

68000 

66250 

.060 

70000 

68272 

.054 

72000 

70292 

.009 

74000 

72312 
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Table  3.  Range  sector  probadsility  of  detection  estiaateo  for  a 
lambda-sigma  jimp  process,  a  lateral  range  of  10,000  meters,  a 
relaxation  time  t  of  3  minutes  and  a  3-out-of-5  criterion. 
Rang©  sector  and  distance  from  CPA  values  are  in  meters  and  zero 
probability  values  are  not  listed. 


prolsabllity 

range  sector 

distan 

.01 

76000 

-74330 

.147 

74000 

-72312 

.090 

72000 

-70292 

.061 

70000 

-68272 

.046 

68000 

-66250 

.033 

66000 

-64226 

.017 

64000 

-62201 

.005 

62000 

-60175 

.001 

60000 

-58146 

.003 

52000 

-50010 

.002 

50000 

-47969 

.001 

48000 

-45924 

.003 

46000 

-43875 

.002 

42000 

-39762 

.001 

40000 

-37697 

.001 

38000 

-35623 

.001 

28000 

-25080 

.001 

26000 

-22913 

.001 

24000 

-20712 

.003 

20000 

-16155 

.002 

18000 

-13748 

.012 

12000 

-4583 

.008 

10000 

0 

.001 

14000 

8307 

.001 

16000 

11180 

.003 

20000 

16155 

.001 

22000 

18466 

.001 

26000 

22913 

.003 

28000 

25080 

.002 

30000 

27221 

.001 

32000 

29343 

.001 

36000 

35541 

.001 

40000 

37696 

.001 

44000 

41821 

.003 

46000 

43875 

.004 

48000 

45924 

.002 

50000 

47969 

.001 

56000 

5408J 

.005 

5BOOO 

56116 

.005 

60000 

58146 

.010 

62000 

60175 

.040 

64000 

62201 

.029 

66000 

64226 

.039 

68000 

66250 

.04j 

70000 

68272 

.046 

72000 

70292 

.005 

74000 

72312 

25 


Table  4.  Range  sector  conditional  probability  of  detection 
estisiates  for  a  gauss-iaarkov  process,  a  lateral  r^ge  of  10,000 
meters,  a  relaxation  time  of  3  minutes  and  a  3-out-of-5 
criterion.  Mnge  sector  and  distzmce  from  CPA  values  are  in 
meters  and  zero  probability  values  are  not  listed. 


probability 

rcmge  sector 

dlstani 

.01 

76000 

-74330 

,133 

74000 

-72312 

.104 

72000 

-70292 

.081 

70000 

-68272 

.041 

68000 

-66250 

.066 

66000 

-64226 

.011 

64000 

-62201 

.005 

62000 

-60175 

.002 

50000 

-47969 

.002 

48000 

-45924 

,002 

40000 

-37697 

.003 

38000 

-35623 

.002 

26000 

-22913 

.011 

12000 

-4583 

.011 

10000 

0 

.002 

12000 

4583 

.002 

32000 

29343 

.002 

38000 

35623 

.002 

46000 

43875 

.003 

50000 

47969 

.002 

52000 

50010 

.002 

58000 

56116 

.003 

60000 

58146 

.012 

62000 

60175 

.068 

64000 

62201 

.055 

66000 

64226 

.088 

68000 

66250 

.128 

70000 

68272 

.132 

72000 

*^0292 

.025 

74000 

72312 
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Table  5.  Range  sector  conditional  probability  of  detection 
estimates  for  a  lambda-sigma  jump  process,  a  lateral  range  of 
10,000  meters,  a  relaseation  time  t  of  3  minutes  and  a 
3-out-of-5  criterion.  Range  sector  and  distance  from  CPA  values 
are  in  meters  and  sero  probability  values  are  not  listed. 
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di4itan<^  from  CPA 


prOltability  rangm  sector 


.01 

.148 

.107 

.Oil 

.068 

.051 

.027 

.008 

.001 

.005 

.003 

.002 

.005 

.003 

.002 

.002 

.002 

.002 

.002 

.003 

.005 

c021 

.014 

.002 

.002 

.003 

.002 

.002 

.006 

.004 

.002 

.002 

.002 

.002 

.006 

.003 

.004 

.002 

.010 

.010 

.020 

.079 

.063 

.090 

.109 

.131 

.016 


78000 

74000 

72000 

70000 

88000 

68060 

84000 

82060 

60000 

52000 

50000 

48000 

48000 

42000 

40000 

38000 

28000 

28000 

24000 

20000 

18000 

12000 

10000 

1400C 

16000 

20000 

22000 

26000 

20000 

30000 

32000 

36000 

40000 

44000 

46000 

40000 

50000 

S60Q0 

58000 

60000 

62000 

64000 

66000 

66000 

70000 

72000 

74000 


-743^0 

-72312 

-70292 

-68272 

-66250 

-84226 

-62201 

-60175 

-56146 

-50010 

-47969 

-45924 

-43875 

-39762 

-37697 

-35623 

-25080 

-22913 

-20712 

-16155 

-13748 

-4583 

0 

8307 

11180 

16153 

18466 

22913 

25080 

27221 

29343 

35541 

37696 

41821 

43875 

45924 

47969 

54063 

56116 

58146 

60175 

62201 

64226 

66250 

68272 

70292 

72312 
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Vex  ttM  siaulatcd  •ncountars,  tho  lateral  range  function 
plots  in  Section  ZV  suggest  that  the  probability  of  detection 
during  an  encounter  im  a  saximun  for  a  relaxation  time  in  the 
neighborhood  of  3  ainutes  for  the  gauas-sarkov  process  ac  well  as 
for  the  laad9da-«igaa  jusp  process.  A  heuristic  argusent  for  this 
apparent  relation  between  probability  of  detection,  relzueation 
tine  and  detection  criterion  for  the  laobda-sigma  junp  process  is 
given  in  i^pendix  3. 

For  the  simulated  encounters  and  a  relaxation  time  of  3 
minutes,  the  lateral  range  function  plots  in  Section  IV  also 
suggest  that  the  probability  of  detection  during  an  encounter  for 
the  gauss~aarko/* process  and  for  laabda-sigma  jua^  process  are 
essentially  equivalent.  This  apparent  equivalenco  does  not  hold 
generally.  For  example,  for  relaxation  tiiaas  of  10  minutes  and 
60  minutes,  the  probability  of  detection  during  an  encounter  for 
the  9auss<'aarKov  process  is  greater  than  that  for  the  lambda- 
sigma  jump  process.  However,  for  both  processes,  the  lateral 
range  function  plots  for  a  relaxation  time  of  zero  correspond  to 
encounters  in  %diich  the  random  ec»aponent£  of  the  signal  excess 
are  independent  gaussian  random  variables.  And,  for  both 
processes,  the  lateral  range  function  plots  for  a  relaxation  time 
of  infinity  correspond  tc.  encounters  in  which  the  random 
components  of  the  signal  excess  are  dependent  gaussian  remdom 
variables  with  correlation  coefficients  of  one.  Consequently,  in 
the  limit  for  a  relaxation  tine  of  zero  or  infinity,  the  two 


ptoeiasM  ar«  a^valant.  Flgura  5  ahcwe  their  lateral  range 
pieta  fdr  a  relaxation  tliui  of  aero  ai^  Figure  6  shows  their 
lateral  range  plota  for  a  relaxation  tine  of  infinity. 

In  ¥able  6  on  Fage  33,  for  the  noeinal  encounter  ping  cycle 
tiae  of  60  aeconda  and  a  laohda-oigaa  juap  process,  conditional 
prd&ahility  of  detection  estiBatea  are  listed  for  a  l-out>o£-l 
detection  criterion  and  a  3-cut-of-5  detection  criterion,  (in 
a  1-out-of-l  detection  criterion,  a  return  is  unconditionally 
classified  ae  an  echo  if  and  only  if  the  signal  excess  associated 
with  the  return  is  greater  than  or  equal  to  eero.) 

Also,  in  Table  6,  conditional  probabilities  of  detection  are 
listed  for  an  encounter  with  a  single  ping  per  sector,  l-out-of>l 
detection  criterion.  For  an  encounter  with  a  single  ping  per 
sector,  l-out-of-1  detection  criterion,  f(SB/o)  is  the 
coiuiitionaX  probability  of  detection  for  a  range  sector  where  « 
indicates  the  standard  noraal  (gauasian)  cunulatlve  distribution 
function.  Conditional  probabilities  of  detection  deteminad  for 
a  range  sector  by  4(SE/a}  have  been  used  to  construct 
probability  of  detection  aaps.  Table  6  suggests  that  for  a 
lasbda-eigaa  jiaap  process,  the  conditional  probabilities  of 
detection  for  a  3-out-of'-5  detection  criterion  nay  not  be 
significantly  different  froa  the  values  detenained  by  ♦(SE/o). 
However,  in  Appendix  3,  a  heuristic  argument  indicates  that  for  a 
laabda-sigaa  juap  process,  a  3-out-of-5  encounter  conditional 
probability  of  detection  estioate  is  of  the  order  of 
#(SE/<J)  •  (t/T)  •exp(-8/T)  where  s  is  the  time  required  for  3 


o^Eis«Otttlv«  succossM  to  occuT  ond  for  this  auii>  t  o&n  be 
cofiisidered  to  be  the  time  in  a  range  sector.  Fox  the  encounter 
elaulation  that  generated  the  data  in  i’able  6,  both  s  and  t 
are  egual  to  3  ainutea  and|  for  the  range  oectors  with  reasonable 
saaple  sizes,  t/v  is  a^rorinately  equal  to  3.  This  iiaplies 
that  the  3>ottt*-of-5  values  should  be  of  the  order  of  «(SE/a) 
idiich  is  consistent  with  Table  6.  For  £in  encounter  whose  lateral 
range  is  10  kilc^eters,  the  length  of  the  range  sector  track 
segnents  go  fron  zero  for  a  range  of  lOO  kilometers  to  l32vS6 
liters  for  the  CPk  range  of  10  kilometers.  For  the  nominal  ping 
cycle  time  of  60  eeconds,  the  number  of  pings  in  the  cpjk  rang  a 
sector  track  segment  ie  approicimately  42.  A  table  similar  to 
Table  6  is  generated  for  the  gauss-markov  process  for  the  seune 
encounter  conditions. 

Table  6  also  suggests  that,  for  a  three  minute  relaxation 
time,  the  1-out-of-l  criterion  is  juperior  to  the  3-out-of-5 
criterion  because  of  the  increased  prt^^ility  of  detection  and 
the  minimum  decision  delay.  However,  for  an  operational  sonar 
system,  tlie  false  alarm  probability  would  be  significantly  lower 
for  the  3-out-of-5  criterion  given  typical  single  return  false 
alarm  probabilities. 
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cingle  ping 

1-out-cf-l 

3-out-of-5 

range  sector 

.000 

76000 

.017 

.047 

.015 

74000 

.130  * 

.323 

.127 

72000 

.130 

.192 

.107 

70000 

.077 

.180 

.079 

63000 

.070 

.119 

.055 

66000 

.067 

.149 

.064 

64000 

.032 

.028 

.022 

62000 

.008 

.026 

.010 

60000 

.002 

.003 

.002 

58000 

.001 

.007 

.002 

56000 

.001 

.007 

54000 

.002 

.007 

.002 

52000 

.004 

.017 

.005 

50000 

.005 

.003 

.005 

48000 

.005 

.017 

.003 

46000 

.003 

.003 

44000 

.002 

.007 

42000 

.002 

40000 

.00? 

.002 

38000 

.002 

.007 

36000 

.002 

.011 

.002 

34000 

.002 

.003 

.002 

32000 

.001 

30000 

.002 

.007 

28000 

.002 

.011 

.002 

26000 

.002 

.011 

24000 

.001 

22000 

,001 

.004 

20000 

.001 

.005 

18000 

.002 

.007 

16000 

.002 

.007 

14000 

.002 

.008 

.005 

12000 

.007 

.111 

.022 

10000 

Table  6,  Range  sector  conditional  probabilities  of  detection  for 
a  single  ping  per  sector,  1-out-of-l  criterion  and  conditional 
probability  of  detection  estimates  for  a  lambda-sigma  jump 
process  and  a  1-out-of-l  criterion,  a  3-out-of-5  criterion 
and  a  10,000  meter  lateral  range.  Range  sector  values  are  in 
meteirs  and  zero  probability  values  are  not  listed. 
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A|^jp«udlx  1.  Program  OaiKuriptione 

The  programs  that  vere  used  to  compute  encounter  detection 
prcdsahility  estimates  are  listed  in  Appendix  2.  The  programs  are 
written  in  BASICA  and  vere  run  on  an  IBK-PC  using  QuickBasic 
which  is  licensed  by  Microsoft  Corporation. 

When  a  program  is  run,  em  initial  prompt  provides  a  choice 
of  either  generating  a  file  containing  encounter  probability 
estimates  or  printing  the  data  contained  in  a  previously 
generated  file.  When  the  progrem  is  used  to  generate  encounter 
probability  estimates,  additional  prompts  regueat  parameter 
values  that  eure  required  to  define  an  encounter. 


2.  Progreua  Listings 

10  SSM  active  gauss-ru&rkov  lateral  rtmge  function/range  sector 
proisabilities 
20 

30  DIM  TL(50) ,PP(20) ,N0(20) ,LR(20) ,CTA(20) ,M1(150) ,P1(150,20) , 
P2C150,20);D1(I50,20),Y1(150,20) 

40  FX*^*AITN(1) 

50  iKfDT  "gensraue  file  or  print  file  (g/p)";N$ 

60  IF  OR  N$-"g"  THEN  GOTO  230 

70  IF  OP  H$>»«p«  THEN  GOTO  80  ELSE  GOTO  50 

80  PRINT 

90  ON  ERROR  GOTO  100:  GOTO  110 

100  RESUME  90 

110  INPUT  <»file  nasa”;A$ 

120  OPEN  «I",#1,A$ 

130  IMPUTll ,RO, SL,TS , NL, DI , TAU, SI6, W, LRM, ST, CT,N3 , NO , NS ,M8 , B$ , 

P$  L 

140  FOR  J-1  TO  L 

150  INP0T#1,  PP(J),LR(J),M1(J),N0(JJ,CTA(J) 

160  FOR  I«0  TO  K1(J} 

170  IKPUT#1,  P1(I,J),R1(I,J),Y1{I,J) 

180  Nmrr  I 
190  NEXT  J 
200  CLOSE 

210  ON  ERROR  GOTO  0 
220  GOTO  1440 
230  PRINT 

240  INPUT  •’aaxiaim  tranamission  loss  lango  (50000/100000) "  ;HTLR 
230  IF  KTLR-50000!  THEN  RI-1000  ELSE  IP  HTLR-IOOOOO!  THEN  RI=2000 
ELSE  GOTO  230 

260  INPUT  "transmission  loss  data  entry  (K*keyboard/D»»disk) "  ?A$ 

270  IP  A$«»"D"  OR  A$»"a"  THEN  GOTO  400 

280  IF  A$»"K«  OR  A$»"k"  THEN  GOTO  290  ELSE  GOTO  260 

290  FOR  I»1  TO  50 

300  A$*i^TR$  (1*2000) 

310  PRINT  ^transmission  loss  at  "A$"  meters  (dB)":  INPUT  TL(I) 

320  JJEST  I 
330  TL(0)-TL(1) 

340  INPUT  "transmission  loss  file  naae";B$ 

350  OPEN  "0«J1,,B$ 

360  FOR  I«0  TO  iO 
370  NRITEfl,  TL(I) 

380  NEXT  I 

390  CLOSE;  GOTO  490 

400  ON  BRROl?  GOTO  410:  GOTO  420 

410  RESUME  <00 

420  INPUT  "transmission  loss  file  naHG";B$ 

430  OPEN 

440  FOR  I»0  TO  50 
450  INPU*;’#1,  TL(I) 

460  NEXT  I 
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470  GSoOSB 

480  OH  EHKOR  GOTO  0 

490  XHHTT  "recognition  differential  (dB)'';HD 

500  msm:  "source  level  (dB)";SL:  INPOT  "target  strength  mB)";TS 
510  IHPOT  "noise  level  (dB)";NL:  INPUT  "directivity  index  (dB)"; 

520  FCHC^SIH^S-NIH’DI-RD:  R£K  figure  of  aerit 

530  XHPOT  "relaxation  tiae  tau  (ainutec,  -1  for  infinity)  "?TAU 

540  INPUT  "sigaa  (dB)"?SIG 

550  INPUT  "relative  spaed  (Jcnots)";?? 

560  IP  N<“0  TESN  PRINT  "Bust  be  greater  than  0";  GOTO  550 
570  H<47<>18S2/69:  RSI  relative  speed  in  Beters  per  ainute 
980  INPUT  "oaxiBUB  lateral  r2uige  (Beters)  "rUtK 
590  IP  LRi0.8"HTIiR  THEN  PRIivi*  "aaxiauB  is  * 

+STR$(.8*}{TIJl)+"  neters":  GOTO  580 

600  INPOT  "lateral  range  step  (neters) "j ST 

610  IP  ST>LRM  THEN  PRINT  "naxiauB  step  is  "+STR$(LRM):  GOTO  600 
620  I^INT(liM/ST) 

630  IP  L>20  THEN  PRINT  "BiniBUB  step  is  "+STR${I^!V20)  :  GOTO  600 

640  INPUT  "ping  cycle  tiae  (seconds) ‘^;CT:  CTM-300 

650  IP  CT>CTH  THEN  PRINT  "Bust  be  less  than  “+STR$(CTM)+ 

"  seconds":  GOTO  640 

660  CT»=CT/60:  REH  ping  cycle  tiae  in  ainutes 

670  N5-INT(SQR(HTLR*IfrLR~LH}»*JJlM)/W/CT) :  REH  number  of  pings  to 
CPA 

680  INPUT  "pings  in  detection  criterion  window" ?N8 

690  IP  Na>*2*N5  THEN  PRINT  "must  be  less  than  "+STR$(2"N5) : 

GOTO  680 

700  IF  N8<1  THEN  PRINT  "ainimuB  is  1  :  GOTO  480 

710  INPUT  "echoes  required  for  detection" ?M8 

720  IP  H8<1  THEN  PRINT  "minimum  is  1":  GOTO  710 

730  IP  M8>N8  THEN  PRINT  "maximua  is  "+STR$(N8)t  GOTO  710 

740  INPUT  "repetitions"  ;K3 

750  IF  N3<“0  THEN  PRINT  "must  be  greater  than  0"t  GOTO  740 
760  INPUT  "data  file  name";A$ 

770  LR-0 

780  FOR  J«1  TO  L 

790  LR*LR-fST:  LR(J)"*LR:  REH  lateral  range 

800  ¥0-SQR(HTLR*KTLR-LR*LR) :  HIM  distance  from  CPA  at  encounter 
start 

810  TE**Y0/W:  REM  TB*tia©  to  CPA  in  minutes 

820  N0»INT{TE/CT) :  IF  N0*0  THEN  N0»1;  REK  adjusted  number  of 

pings  to  CPA 

830  CTA*TE/N0:  REH  adjusted  ping  cycle  tiES  in  ainutes 
840  DIM  SE(2"NO),P(2*NO},Y{2*NO),R(2*NO),N(2*NO) 

850  TK— TE 

060  FOR  K»0  TO  NO 

870  Ri-SQR(LS*tR+W*W*TK*TK) :  IK°INT(R1/RI)  REM  transmission  loss 

range  index 

880  SE«FOH-2*TL(rN) 

890  SE(K)°SE;  TK-*TK-fCTA:  R(K)=IN*RI  REH  transmission  loss  ranae 
for  POM 


900  NEXT  K 

910  FOR  K*N0+1  TO  2*N0 

920  SE(K)»SE(2*N0~K):  R(K)«R(2*N0“K) 

930  NEXT  K 
940  P«0 

950  FOR  I»1  TO  »3 
960  DIM  M(N8“1) 

970  OOSUB  2120 
980  X»SIG*Q 

990  IF  TAU  *“1  a>rO  1020 

1000  IF  TAU®0  THEN  RHO=0:  FAC=1:  GOTO  1020 

1010  RHO«EXP(»CTA/TAU) :  FAC=SQR{l-RHO*RHO) ;  REM  gauss-jaarkov 

parameters 

1020  FOR  K=0  TO  2*N0 

1030  XF  K*0  OR  TAU^-l  GCTO  1060 

1040  GOSUB  2120 

1050  X»SIG*FAC*Q+RHO*X 

1060  XSE=SE(K)+X 

1070  IF  K<*N8-1  THE!?  J1=Kj  GOTO  1120 
1080  FOR  J1=0  TO  N8-2 
1090  M(J1)«*M(J1+1) 

1100  NEXT  J1 
1110  J1-N8-1 

1120  IF  XSE>=0  THEN  M(J1)»1  ELSE  M(J1)=0 
1130  S«0 

1140  FOR  J1=0  TO  N8-1 
1150  S-S+M(J1) 

1160  NEXT  J1 

1170  IF  S>«M8  GOTO  1200 

1180  NEXT  K 

1190  GOTO  1210 

1200  P-P+1:  P(K)»P(K)+1 

1210  SR2^E  M 

1220  NEXT  I 

1230  PP(J)=P/N3 

1240  P1(0,J)=P(0) :  R1(0,J)=R(0) ;  Y1(0,J)=-Y0 
1250  FOR  I«1  TO  2*N0 

1260  IF  R(I)=R(I-1)  THEN  N(I)=N(I-1) :  P1(N(I) , J)=P(I)+P1(N(I) , J) 
GOTO  1310 

1270  N(I)>=N(I-l)+l:  P1(N(I)  ,J)=P(I)  :  R1  (N(I)  ,  J)  =R(I) 

1280  IF  LR>R1(N(I) ,J)  THEN  Y1(N(I) , J)=0;  GOTO  1310 
1290  Y1(N(I) ,J)=SQR(R1(N(I) ,J)*R1(N(I) ,J)-LR*LR) 

1300  IF  Y1(N(I) ,J)<ABS(Y1(N(I-1) ,J))  THEN  Y1(N(I) , J)=-Y1(N(I) , J) 
1310  NEXT  I 

1320  M1(J}“N(2*N0) :  N0(J)=N0;  CTA(J) =CTA*60 

1330  ERASE  P,Y,R,N,SE 

1340  NEXT  J 

1350  OPEN  ''0«,#1,A$ 

1360  NEITE#1,RD,SL,TS,NL,DI,TAU,SIG,W,LRM,ST,CT,N3,N0,N8,M8, 
B$,P$,L 

1370  FOR  J«1  TO  L 

1380  WEITEn,  ?P(J)  ,M1(J)  ,N0(J)  ,CTA(j) 


\p$ 

%B$ 

N 


1390  FOR  1*0  TO  M1(J) 

1400  WRITE#1,  P1(I,J) ,R1(I,J) 

1410  NEXT  I 
1420  NEXT  J 
1430  CLOSE 
1440  PRINT 

1450  INPUT  "print  encounter  pareuaeters  (Y*yes/N-no) "?D$ 

1460  IP  D$-"Y"  OR  D$-"y«  GOTO  1480 
1470  IP  D$«"N"  OR  D$-"n"  GOTO  1660  ELSE  GOTO  1450 
1480  LPRINT 

1490  LPRINT  "program  file  name 
1500  LPRINT  "transmission  loss  file  name 
1510  LPRINT  "recognition  differential  (dB) 

1520  LPRINT  "source  level  (d3) 

1530  LPRINT  "target  strengtt  (dB) 

1540  LPRINT  "noise  level  (dB) 

1550  LPRINT  "directivity  index  (dB) 

1560  LPRINT  "relaxation  time  tau  (minutes) 

1570  LPRINT  "sigma  (dB) 

1580  LPRINT  "relative  speed  (knots) 

1590  LPRINT  "xaaximum  lateral  range  (meters) 

1600  LPRINT  “lateral  rauige  step  (meters) 

1610  LPRINT  "ping  cycle  time  (seconds) 

1620  LPRINT  "pings  in  detection  criterion  window 
1630  LPRINT  "echoes  required  for  detection 
1640  LPRINT  "repetitions 
1650  LPRINT  "data  file  name 
1660  PRINT 

1670  INPUT  "print  encounter  probabilities  (Y»»yes/N«no) "  ;D$ 

1680  IF  D$-"Y"  OR  D$-"y"  GOTO  1700 
1690  IF  D$“"N"  OR  D$»"n"  GOTO  1780  ELSE  GOTO  1670 
1700  LPRINT:  LPRINT:  LPRINT:  LPRINT  A$+"  encounter  probabilities" 
1710  FOR  J>«1  TO  L 
1720  LPRINT:  LPRINT 
1730  LPRINT  "lateral  range 
1740  LPRINT  "probability  of  detection 
1750  TaPRINT  "adjusted  ping  cycle  time  (seconds) 

1760  LPr.INT 
1770  NEXT  J 
1780  PRINT 

1790  INPUT  "print  range  sector  probabilities  (Y=yes/N»no) " ;D$ 

1800  IF  D$«"Y"  OR  D$»«y«  GOTO  1820 
1810  IF  D$“"N"  OR  GOTO  1930  ELSE  GOTO  1790 

1820  LPRINT:  LPRINT:  LPRINT:  LPRINT  A$+«  range  sector 
probabilities" 

1830  FOR  J=1  TO  L 
184-0  LPRIirr;  LPRINT 
1850  LPRINT  "lateral  range 
1860  LPRINT  "probability  of  detection 
1S70  LPRINT  "adjusted  ping  cycle  time  (seconds) 

1880  LPRINT 
1890  FOR  I«0  TO  M1(J) 


,RD 
",SL 
",TS 
",NL 
",DI 
",TAU 
",SIG 

",W"60/1852 

",LRM 

",ST 

«,CT*60 

",N8 

",M8 

",N3 

",A$ 


M«(J) 

",PP(J) 

",CTA(J) 


MmR(J) 

",PP(J) 

",CTA(J) 


1900  LPRIH?  “PC"!")  -  *'P1(I,J)/N3  TAB(30)  -  "R1(I,J) 

TAB(60)  "YC"!")  -  ••Y1(I,J) 

1910  NEXT  I 
1920  NEXT  J 
1930  PRINT 

1940  INPUT  "print  range  sector  conditional  probabilities 
(Y*yes/N»^o) » yD$ 

1950  IF  D$-"Y"  OR  D$-"y"  GOTO  1970 

I960  IP  D$»"N"  OR  D$-"n"  GOTO  2110  ELSE  GOTO  1940 

1970  LPRINT:  LPRINT:  LPRIFT:  LPRINT  A$+"  range  sector  conditional 

probabilities" 

1980  FOR  J-1  TO  L 
1990  IjPRINT:  LPRINT 

2000  I^RINT  "lateral  range  ",LR(J) 

2010  LPRINT  "probability  of  detection  ",PP(J) 

2020  LPRINT  "adjusted  ping  cycle  time  (seconds)  ",CTA(J) 

2030  LPRINT 
2040  D-N3 

2050  FOR  I-O  TO  M1(J) 

2060  IF  D-0  THEN  P2(I,J)*0  ELSE  P2 (I, J)»P1(I, J)/D 

2070  LPRINT  "P(«I«)  «  *P2(I,J)  TAB(30)  "R("I")  "R1(I,J)  TAB(60) 

"Y(«I")  -  "Y1(I,J) 

2080  IF  D>0  THEN  D-D-P1(I,J) 

2090  NEXT  I 
2100  NEXT  J 
2110  END 

2120  RV«PJ?D;  IF  RV-0  THEN  RV’2.8E“38 
2130  IF  RV-1  THIN  RV“l-2.8E-38 

2140  C2«SQR(-2*LOG(RV))*SIN(2*PI«KND) :  REH  standard  normal  random 

number 

2150  RETURN 


10  B8N  60tlv«  le»bda~sigaa  lateral  rangfa  funotlon/ranga  aact.or 

|>rO£}abilitias 

20  P$»*'RSLSA.Bi^** 

30  DIM  TL(SO)  ,PF(20)  ,N0(20)  ,ZA(20)  ,CTA(2C)  ,!£1(150)  ,P1(150,20) , 

P2 (150,20) ,R1 (150, 20) ,¥1(150, 20) 

40  PI-4*ATH(1):  RCla.88623:  RC2-. 21459 
50  ZHHJT  "ganarata  flla  or  print  flla  (g/p)''?N$ 

60  IP  N$«"G"  OR  N$»*g''  THEN  GOTO  230 

70  IF  OR  THEN  GOTO  80  ELSE  GOTO  50 

60  PRINT 

90  ON  ERROR  GOTO  100:  GOTO  110 

100  RESOME  90 

110  INPUT  "file  naaa“;A$ 

120  OPEN  ®*I'»,#1;A$ 

130  INPOT#l,  RD,SL,TS,NL,DI,TAD,SIG,W,1RM,ST,CT,N3,NC,N8,MQ, 

B$,C$,P$,L 

140  POR  J»1  TO  L 

150  INPUT#1,  PP(J),LR(J),M1(J),N0(J),CTA(J) 

160  FOR  I«0  TO  K1(J) 

170  INPUT#1,  P1(I,J),R1(1,J),Y1(I,J) 

130  NEXT  I 
190  NEXT  J 
200  CLOSE 

210  ON  ERROR  GOTO  0 
220  GOTO  1520 
230  PRIN^T 

243  IN^JT  "aaxiaim  trananlsalon  loss  range  (50000/100000)  **  ;MTLR 
250  IP  MTI^-SOOOOl  THEN  RI«-1000  ELSE  IP  HTIil-1000001  THEN  Rl-2000 
ELSE  GOTO  230 

260  INPUT  "transHisaioa  loss  data  entry  (K-koyboard/D«dislc)  “?A$ 

270  IP  A^-^O"  OR  A$-*»d“  GOTO  400 

280  IF  A$-"K*  QE  A$«“k»  Gi^SQ  290  ELSE  GOTO  260 

290  POE  I»1  TO  50 

300  A$»£TR$ (1*2000) 

310  PRINT  "transaisssion  loe®  at  •*A$"  aeters  (dB)":  INKJT  TL(I) 

320  NEXT  I 
330  TL(0)»TL(1) 

340  INPUT  “transaiaaion  loss  file  naae";S$ 

350  OPEN  B$ 

36C  FOR  I»0  TO  50 
370  WEITEII,  TL(I) 

380  NEXT  I 

390  CLOSE:  GtOTO  520 

400  OR  ERROR  GOTO  410:  GOTO  420 

410  RESUME  400 

426  INPUI'  *'transjQission  loss  file  nei3Q";BS 
430  OPEN  “I",#!,  BS 
440  FOR  I«0  TO  SO 
450  INPUTfl,  TL(I) 

460  NEXT  I 
470  CLOSE 

480  ON  ERROR  GOTO  0 
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490  "aignal  esecus  distribution  (Gxgaussian/^raylsigh)  ";c$ 

500  IF  C$»«6"  OR  C$«*'g*'  GOTO  520 

510  XF  C$»"R«  OR  C$»**r«  GOTO  520  ELSE  GOTO  490 

520  IMH3T  "racognition  diffarantial  (dB)";R0 

530  XHFCT  *aouroa  laval  (dB)*;SL:  IHFtlT  "target  strength  (dB)";TS 
540  XBFOT  "noise  level  (dB)";ML:  IHFUT  "directivity  index  (dB)"; 
DI 

550  FOif-SLfTS-HIrfDZ-RDs  REM  figure  of  nerit 

560  ZNFOT  "relaxation  tisa  tau  (ainutes,  -1  for  inf inity) "  ;TAtJ 

579  ZMSHT  "signa  (dB)"}SX6 

580  IMFOT  "relative  spaed  (knots) ";W 

590  IF  M<"*0  THEN  PRINT  "must  be  greater  than  0":  y^OTO  580 
600  11^*1852/60;  REK  relative  speed  in  neters  per  minute 
610  2MP0T  "saxiaua  lateral  range  (meters) ";LRM 
620  IF  LRH>.8*MT1II  THEN  PRINT  "maximum  is  "-fSTROC  .8*HTI^)-i- 
"  meters":  GOTO  610 

630  INPUT  "lateral  range  step  (meters)"; ST 

640  IF  ST>IRH  THEM  PRINT  "maximum  step  is  "+STR$(LEH)  :  GOTO  630 
650  I^IMT(IiRN/ST) 

660  IF  L>20  THEM  PRINT  "minimum  step  is  "-»-STR$(IRII/20)'f"  meters": 
GOTO  630 

670  INPUT  "ping  cycle  time  (seconds) ";CT:  CTM»300 
636  IF  CT>CTK  THEN  PRINT  "must  be  less  than  "-fSTRO  (CTN) 

"  eeconda":  GOTO  670 

690  CTK^T/SO:  REH  ping  cycle  time  in  minutes 

700  N5«INT(SQR(MTLR*MTIR-IRM*LR24)/N/cr) ;  REM  number  of  pings  to 
CPA 

710  INPUT  "pings  in  detection  criterion  window" ;N8 

720  IF  N3>«2*N5  THEM  PRI!iT  "must  be  less  than  "+STR$(2*N5)  ; 

GOTO  710 

730  IF  N8<1  THEN  PRINT  "minimu®  is  I":  SOTO  670 

740  INPUT  "echoes  required  for  detection"  ;!!8 

750  IF  M8<1  THEN  PRINT  "minimum  is  1":  GOTO  740 

760  IF  He>N8  THEN  PRINT  “maximum  is  "+S1R$(N8):  GOTO  740 

770  INPUT  "repetitions" ;N3 

780  IF  H3<»0  THEN  PRINT  "must  be  greater  than  0":  GOTO  770 
790  INPUT  "date  file  naae";A$ 

800  LR«0 

810  rm  J*1  TO  L 

820  UC-LR4ST:  LR(J)-LR:  REH  lateral  range 

S30  YC«SQR(HTIsR*KTIiR-LR*LR) :  REM  distance  from  CPA  at  encounter 
start 

840  T£t»V'0/t?:  REK  TE“tiae  to  CPA  in  minutes 

350  N0-INT{TE/Cr) :  IP  HO-0  THEN  NO-1;  REM  adjusted  number  of 
pings  to  CPA 

860  CTA-TE/NO:  RES  adjusted  ping  cycle  tine  in  minutes 
870  DIN  SB(2«K0) ,P(2*N0) ,y(2*N0) ,R(2*N0) ,N{2*N0) 

880  TK— TB 

390  FOR  K-0  TO  NO 

900  Rl-SQR(Iil*LR+W*W*TK*TK) ;  IN-INT(R1/RI)  TL  range  index 
910  RE-POM-2 *TL( IN) 
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9ao  8B(K}«6B:  ¥K««KHri!As  R(X)«ZM*RZ  BBX  tranaaiMlon  Iom  ranga 

for  fm 

030  MttlWP  V 

$40  FOR  TO  2*H0 

9$0  8B(X)»SS(2*H0-K):  R(K)«R(2*H0-K} 

960  B6XT  K 
970  F-0 

980  FOR  Z«I  TO  N3 
990  DIN  M(Ne-l) 

1000  IF  TRU«0  GOTO  1050 
1010  GOSUB  2210 
1020  IF  TAU—1  GOTO  1050 
1030  GOSUB  2270 

1040  K*">TAO*L06(1-RV)/CTA:  REN  nuabar  of  pings  to  tha  first  jium 

1050  FOR  X«0  TO  2*K0 

1060  IF  TAU-0  GOTO  1110 

1070  IF  K<H  OR  TAI>*-1  GOTO  1120 

1080  GOSUB  2270 

1090  CH— TAU*L0G(1-RV)/CTA 

1100  K«N4'D19:  REN  nuabar  of  pings  to  tha  naxt  juap 
1110  GOSUB  2210 
1120  XSE»SB(K)-H3 

1130  IF  R<«H8-1  THEN  J1«K:  GOTO  1180 
1140  FOR  Jl-0  TO  N6-2 
1150  H(J1)«N(J14-1) 

1160  NEXT  J1 
1170  J1-N8-1 

1180  IF  XSB>«0  THEN  H(J1}«1  ELSE  N(J1)»0 
1190  S«0 

1200  FOE  Jl-0  TO  N8-1 
1210  S*S+M(J1) 

1220  HBXT  J1 

1230  IF  S>«iN@  THEE  GOTO  1260 

1240  NEXT  R 

1250  GOTO  1270 

1260  F-F+ls  P(K)«P(K)+1 

1Z70  ERASE  H 

1280  NEXT  I 

1290  PP{J)-P/N3 

1300  K<0)«0:  Pl(0,a)«P(0) s  R1(0,J)»R(0) ;  V1(0,J)«-V0 
1310  FOR  I-l  TO  2*N0 

1320  IF  R(I)-R(I-1)  T^HES  K(I)»R(I-1);  ?1{N(I)  ,J)-P(I)+P1(N(I)  ,J) 
GOTO  1370 

1330  H(I)-N(I-l)+l:  P1{N(I)  ,J)-P{I)  :  R1(N(I) ,  J)«.R(I) 

1340  IF  LR>R1(H(I),J)  THEN  yi(K(Z) , J)-0;  GOTO  1370 
1350  V1(S(I) ,J)»SQR(R1(N(I) ,J)*R1(N(I) ,J)-IE«LR) 

1360  IF  Yl(ll{I),J)<ABS{Yl(N(I-l),J))  THEN  V1{N{I)  , J)-- Y1(H(I)  ,J) 
1370  NEXT  I 

1380  K1CJ)«N(2*H0) :  N0(J)»N0:  CTA(J)»CTA*60 
1390  ERASE  P,V,R^N,SE 
1400  NEST  J 

1410  IF  OR  c$-*g«  THEN  C$-"gaussian" :  GOTO  1430 


1420  C$»'*rayl«l9h'' 

1430  OPEH  *0'',«X,A$ 

1440  WSITBfl,  HD,SL,TS,HL,DX,TAD,SI6,W,LEiM,S^.CT,N3,N0,K8, 

K[8,B$,C$,I»$,L 

1450  FOB  J«1  TO  L 

1460  VRITEfl,  PP(J),XA(J),H1(J),H0(J),CTA(J) 

1470  FOE  I-O  TO  M1(J) 

1480  KIElTSfl,  P1(Z,J),R1(1,J),V1(I,J) 

14S0  liSXT  I 
1500  NEXT  J 
1510  CLOSE 
1520  FRIKT 

1530  ZHPtJT  "print  ancountar  p&raaatara  ( £»yes/ Keno) " ;  D$ 

1540  IF  D$-«Y"  OR  D$«"y"  GOTO  1560 

1550  IF  D$-"N"  OR  GOTO  1750  ELSj?  GOTO  1530 


1560  LFRIKT 

1570  LKtIHT  "progran  fila  nase  ",p$ 

1580  LERINT  "tranialaslon  loss  file  ru>me  '%B$ 

1590  LFRIRT  "signal  excess  distribution  ",C$ 

1600  LKlDiT  "recognition  differential  (dB)  ",RD 

1610  LFRINT  "source  level  (dB)  ",SL 

1620  LinsiKT  "target  strength  (dB)  ",TS 

1630  LXRIHT  "noise  level  (dB)  ",£}L 

1640  LPRIHT  "directivity  index  (dB)  ",DI 

1650  L]^NT  "relaxation  tisa  tav  (sinutes)  ",TAn* 

1660  LPRINT  "sigaa  (dB)  ",SI6 

1670  LERIKT  "relative  spaed  (teiits)  ",l?*60/1852 

1680  LXRIKT  "Kaxisun  lateral  range  (asters)  *,LBH 

1690  L^UNT  "lateral  range  step  (aeT>ers)  ",ST 

1700  LFRIliT  "ping  cycle  time  (seconds)  ",CT*60 

1710  LFEUINT  "pings  in  det«5Ction  criterion  window  ",N8 

1720  L^llirr  "echoes  required  for  detection  ",He 

1730  LIRINT  "repetitions  •*,N3 

1740  LPRIIJT  "data  file  nam  »,A$ 

1750  PSIKT 


1760  IKPtJT  “print  encounter  probabilities  <Y“yes/N«no)  "70$ 

1770  IP  D$-«Y«  OS  D$w"T*  GOTC*  1/90 

1780  IP  D$»"N"  OS  D$-«^  30T0  1870  ELSE  GOTO  1770 

1790  LPSINT!  LFSXNT;  LFRIilTi  LPRIIR;  A$+"  encounter  probabilities" 

1800  FOR  J-1  TJ  L 

1810  LPRINT:  I^RINT 

1820  LPRINT  "lateral  range  ",LR(J) 

1830  LPRINT  "probability  of  detection  ",PP(J) 

1840  I^RINT  "adjusted  ping  cycle  tixi©  (seconds)  *,CTA(J) 

1850  LPRINT 
1860  NEXT  Ci 
1870  PRINT 

1380  INPUT  "print  range  sector  probabilities  (Y-yes/N-no) * ;D$ 

1890  IF  D$»"Y"  OR  S$»"y"  Gi^IO  1910 

1900  IP  D$»"N"  OS  (Sol'd  2020  ELSE  GOTO  1880 

1910  LPRINT:  LS^lINT:  LPRINT:  LPRINT  A$+"  rang©  sector 

probabilities" 
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1920  FOR  J-1  TO  L 
1930  ZimilT:  LHIZHT 

1940  LCRX^  "lateral  ranga  ",IR(J) 

1950  iFRXItT  "probability  of  dataction  "«PP(J) 

1960  XiPRIlif  "adjuatad  ping  cycle  tiaa  (aaconds)  ",CTA(J) 

1970  IPRZIIT 

1980  FOR  Z-O  TO  M1(J) 

1990  LHiINT  "PC"!")  -  "P1(I,J)/H3  TAB(30)  "RCI")  »  "RICI^J) 
TAB(60)  "YCI")  -  "Y1(I,J) 

2000  NEXT  I 
2010  NEXT  J 
2020  PRINT 

2030  INPUT  "print  rai^^e  sector  coiulitlonal  probabilities 
CY-y«*/H*no) "  t  D$ 

2040  IF  D$-"Y''  OR  D$-''y"  GOTO  2060 

2050  IF  D$»"N"  OR  D$*"n"  GOTO  2200  ELSE  GOTO  2030 

2060  LPRINT:  LFRINT:  LPRIHT:  LPRIKT  A$4-'*  range  sector  conditional 

probabilities" 

2070  FOR  J«1  TO  L 
2080  LmNT:  LPRINT 

2090  LPRINT  "lateral  range  ",LR(J) 

2100  LPRINT  "probability  of  detection  ",PP(J) 

2110  LPRINT  "adjusted  ping  cycle  tiae  (seconds)  ",CTA(J) 

2120  LFRINT 
2130  0bN3 

2140  FOR  I-O  TO  M1(J) 

2150  IF  D-O  THEN  P2(I,J)-0  ELSE  P2 (I, J)-P1(I, J)/D 

2160  LPRINT  "P("I")  -  "P2(I,J)  TAB(30)  *R("I")  »  "R1(I,J)  TAB(60) 

"Y("I»)  -  "Y1(^,J) 

2170  IP  D>0  THEN  D-D-P1(I,J) 

2180  NEXT  I 
2190  NEXT  J 
2200  END 
2210  GOSUB  2370 

2220  IP  C$«"G«  OK  C$-"g"  GOTO  2250 

2230  Q*SIG*(RC1-SQR(-L0G(RV)))/SQR(RC2) :  REH  randOQ  signal  excess 
part 

2240  RETURN 

2250  Q«SIG*SQR(-2*LOG(RV))"SIN(2*PI*RHD):  REN  randan  signal 
excess  part 
2260  RETURN 

2270  RV-RNO:  IF  RV«0  THEN  RV-i2.8E-38 
2280  IF  RV-1  THEN  RV-1-2.8B-38 
2290  RETURN 


44 


3.  A  Heuristic  Argument 

The  following  hhuristic  argument  applies  to  the  lambda-sigma 
jufi^)  process.  Suppose  that  SE  is  constant  over  a  time  interval 
of  length  t.  Then,  the  occurrence  of  a  jump  during  the  interval 
that  is  sufficient  for  conditionally  classifying  an  echo  as  an 
echo  (a  success)  is  determined  by  a  Poisson  process  with  rate 
«(SE/o)/t.  In  addition,  suppose  *(SE/a)  is  small  enough  so 
that  the  probability  that  one  success  will  occur  in  the  interval 
is  approximately  •(SE/o) • (t/T) .  (If  this  expression  is  equal  to 
.1,  the  probability  of  one  success  in  the  interval  is  .0905  and 
the  probability  of  two  or  more  successes  in  the  interval  is 
.0047.)  The  probability  that  the  time  to  the  first  jump 
following  a  success  will  be  greater  them  or  equal  to  s  is 
exp(-s/T).  ^is  implies  that  the  probability  of  a  success  in 
the  interval  after  which  the  random  component  of  the  signal 
excess  is  constant  for  a  time  at  least  equal  to  s  is 
approximately  ♦(SE/o) • (t/t) •exp(-s/T) .  This  expression  is  a 
maximum  when  t  s.  Let  s  be  the  time  required  for  k 
consecutive  returns  and  let  t  be  a  time  greater  than  but  of  the 
order  of  the  time  required  for  n  consecutive  returns.  Then, 
given  that  the  signal  excess  at  the  beginning  of  a  time  interval 
of  length  t  is  less  than  sero,  the  above  suggests  that  for  a 
k-out-of-n  criterion  the  conditional  probability  of  detection 
during  the  interval  will  have  a  maximum  for  y  in  the 
neighborhood  of  s  if  »(SE/o) • (t/t) •exp(-6/T)  is  small  and  s 
is  sufficiently  smaller  than  t.  This  suggests,  for  t  in  the 
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nelgliborhood  of  s,  that  tha  ranga  aaotor  conditional 
probabilitiaa  for  an  anoountar  will  hava  a  naxtanw  and 
conaaquantly  tha  datactlon  probability  for  tha  ancounter  will 
hava  a  Baxlaua  also. 

As  t  bacomas  saall,  ♦(SE/o) • (t/T) •a3^(-s/T)  will  become 
greater  than  one  and  tha  as^roxlmation  dascrlbad  above  clearly 
does  not  hold.  In  this  case,  giver,  the  signal  excess  at  the 
beginning  of  a  time  interval  of  length  t  is  less  than  zero,  the 
conditional  probability  of  detection  during  an  interval 
approaches  the  probability  of  3  or  more  successes  of  out  of  5 
independent  trials  where  the  probability  of  success  is  fi(SE/a). 

As  T  becomes  large,  the  probability  of  the  random 
component  of  tha  signal  excess  changing  during  an  encounter 
becomes  small  and  in  the  limit  the  random  component  of  the  signal 
excess  will  be  a  constant  value  during  an  encounter  and  the 
signal  excess  can  be  represented  by  SE  +  X  whore  the  value  of 
X  is  determined  at  the  beginning  of  the  encounter.  In  this 
case,  if  SE  +  X  is  greater  than  zero  for  3  returns  out  of  5 
during  an  encounter,  detection  will  occur,  otherwise  it  will  not. 

For  the  encounters  that  were  considered,  for  both  small  and 
large  values  of  t,  the  conditional  prQb^a^ilities  of  detection 
for  the  range  segments  of  an  encounter  appear  to  be  dominated  by 
those  for  T  in  the  neighborhood  of  3  minutes. 
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